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Abstract 

Using the latest determinations of several theoretical and experimental param- 
eters, we update the Unitarity Triangle analysis in the Standard Model. The basic 
experimental constraints come from the measurements of |^fe/V^fe|, Am^, the lower 
limit on Arris, ^K, and the measurement of the phase of the B^-Bd mixing am- 
plitude through the time-dependent CP asymmetry in B^ — > J/t/jK^ decays. In 
addition, we consider the direct determination of a, 7, 2/3+7 cos 2/3 from the 
measurements of new CP-violating quantities, recently performed at the B facto- 
ries. We also discuss the opportunities offered by improving the precision of the 
various physical quantities entering in the determination of the Unitarity Triangle 
parameters. The results and the plots presented in this paper can also be found at 
the URL http://www.utfit.org, where they are continuously updated with the 
newest experimental and theoretical results. 



1 Introduction 



The Standard Model (SM) of electroweak and strong interactions provides an excellent 
description of all observed phenomena in particle physics up to the energies presently 
explored. 

The LEP/SLD era of precision electroweak physics has ended, leaving us a beautiful 
legacy of measurements that strongly constrain alternative mechanisms of electroweak 
symmetry breaking [1], pointing strongly to a light Higgs boson. On the other hand, B 
factories have thoroughly developed the study of precision B physics, which, after the 
great achievement of the Acp{B^ J/ijjK^) measurement, is entering its mature age 
with many analyses aimed at measuring the angles of the Unitarity Triangle (UT) in 
different processes. 

This remarkable experimental progress has been paralleled by many theoretical nov- 
elties: on the one hand, the consolidation of the Heavy Quark Expansion and Lattice 
QCD (LQCD) results on heavy mesons; on the other, the constant refinement of ideas 
and techniques to extract information on the UT angles from B decay rates and CP 
asymmetries. 

Waiting for LHC to uncover the mysteries connected to electroweak symmetry breaking 
and the hierarchy problem, we have the chance to indirectly probe scales up to a few 
TeV through a precision study of flavour physics. The aim of the 'UTfit Collaboration 
is twofold: i) the "classic" issue of assessing our present knowledge of flavour physics in 
the SM, combining all available experimental and theoretical information; ii) the future- 
oriented task of constraining new physics models and their parameters on the ground of 
flavour and CP-violating phenomena. 

Up to now, the standard UT analysis [2, 3, 4] relies on the following measurements: 
iKtfel/lKbl, Amrf, the limit on Am^, and the measurements of CP-violating quantities 
in the kaon [sk) and in the B (sin 2/3) sectors. Inputs to this analysis consist of a 
large body of both experimental measurements and theoretically determined parameters, 
where LQCD calculations play a central role. A careful choice (and a continuous update) 
of the values of these parameters is a prerequisite in this study. The values and errors 
attributed to these parameters in the present study are summarized in Table 1 (Section 2). 
The results of the analysis and the determination of the UT parameters are presented and 
discussed in Section 3 which is an update of similar analyses performed in [2, 3] to which 
the reader can refer for more details. 

New CP-violating quantities have been recently measured by the B factories, allowing 
for the determination of several combinations of UT angles. The measurements of a 
(using TTTT, pp and up modes), 7 (using modes), 2/3-1-7 (using D'^*\[p) modes), 

and cos 2/3 from B^ — > J/ipK* are now available. These measurements and their effect on 
the UT fit are discussed in Section 4. 

Finally in Section 5 we discuss the perspectives opened by improving the precision in 
the measurements of various physical quantities entering the UT analysis. In particular, 
we investigate to which extent future and improved determinations of the experimental 
constraints, such as sin 2/3, Arris, ot and 7, could allow us to invalidate the SM, thus 
signalling the presence of new physics effects. 
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2 Inputs used for the "standard" analysis 

The values and errors of the relevant quantities used in this paper for the standard anal- 
ysis of the CKM parameters (corresponding to the constraints from iKifel/lKb], Am^, 
Arris/ Amd, and sin 2/5) are summarized in Table 1. 

The novelties here are the final LEP/SLD likelihood for Am^, the use of \ Vub\ measure- 
ments from inclusive semileptonic decays at the B factories [5], the updated value of sin 2(5 
and a new treatment of the non-perturbative QCD parameters as explained in the follow- 
ing Section 2.1. In addition, we use updated values of the top mass = 178.0 ± 4.3 
GeV [6] and of the CKM parameter A. The latter comes from the average of the following 
values [7] 

A(K. fromK^a) = 0.2250 ±0.0021 
A(K<i+ unitarity) = 0.2265 ± 0.0020. (1) 

Finally, we now calculate QCD corrections to AB — 2 and AS — 2 processes, fully taking 
into account their correlations with the input parameters and, by varying the matching 
scale, the uncertainty introduced by the residual scale dependence. 



2.1 Use of ^, fs^y Bb^ and fBdyBsd in Am^ and Am^ constraints 

One of the important differences with respect to previous studies is in the use of the 
information from non-perturbative QCD parameters entering the expressions of ArUs and 
Amd- The Bg-Bg mass difference is proportional to the square of the element \Vts\- Up 
to Cabibbo- suppressed corrections, \ Vts\ is independent of p and rj. As a consequence, the 
measurement of Am^ would provide a strong constraint on the non-perturbative QCD 
parameter Jb^Bb^- 

We propose a new and more appropriate way of treating the constraints coming from 
the measurements of Am^ and Am^. In previous analyses, these constraints were imple- 
mented using the following equations 



Am, oc [{l-pr + vVlBsd (2) 
Am, « flM=flaBBd^e 



where ^ = fs^y BbJ fB^y Bb^- In this case, the input quantities are fB^y Bb^ and ^. 
The constraints from Am, and the knowledge of ^ are used to improve the knowledge on 

fsd \J which thus makes the constraint on Am, more effective. 

The hadronic parameter that is better determined from lattice calculations, however, 

is Jb^Bbs, whereas and Bb^ axe affected by larger uncertainties coming from the 
chiral extrapolations. These uncertainties are strongly correlated. For this reason, a 
better approach consists in replacing Eq. (2) with 

Am, « [(l-pf+f]l^ (3) 
Am, oc f^BB,- 
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Parameter 


Value 


Gaussian (a) Uniform 

(half-width) 


A 


0.2258 


0.0014 


\Vh\(exc\ ) 


41.4 X 10~^ 


2.1 X 10"^ 


1 CO 1 


41.6 X 10~^ 


7 X 10"^ 6 X 10^3 

■ 1 ''N -1-\J w « \J XX _L yj 


\V h\(exc\ ) 

1 •'liO 1 ^^C-^V-J-' J 


?,'?, Ox in~^ 


2 4 X 10"'' 4 6 X 10""' 




47.0 X 10~^ 


4.4 X 10~* 




502 ns~^ 


n OOfi n<?~^ 




> 14 5 ns~^ at Q5% C 


oCllol Ul V 1 u V J. O .O 




276 MeV 


38 MeV 




1.24 


0.04 0.06 


Bk 


0.86 


0.06 0.14 


sk 


2.280 X 10-3 


0.013 X 10-3 - 


Ik 


0.159 GeV 


fixc^d 


ISmK 


0.5301 xlO-2 ps-i 


fixed 


sin 2/? 


0.726 


0.037 


mt 


168.5 GeV 


4.1 GeV 


rub 


4.21 GeV 


0.08 GeV 


rric 


1.3 GeV 


0.1 GeV 


as{Mz) 


0.119 


0.003 


Gf 


1.16639 xlO-^GeV-^ 


fixed 


mw 


80.425 GeV 


fixed 


TUbO 

a 


5.279 GeV 


fixed 


TUbo 


5.375 GeV 


fixed 


mi 


0.497648 GeV 


fixed 



Table 1: Values of the relevant quantities used in the UT fit. The Gaussian and the flat 
contributions to the uncertainty are given in the third and fourth columns respectively (for 
details on the statistical treatment see [2]). 



4 



At present, this new parameterization does not have a large effect on the final results. 
It allows, however, to take into account more accurately the uncertainty from the chiral 
extrapolation in lattice calculations of fs^- Note that in this way, in order to obtain a 
more effective constraint on Am^, the error on ^ should also be improved. 

3 Determination of the Unitarity Triangle parame- 
ters: the standard analysis 

Using the constraints from \Vub\ / Am^i, Arris, and sin 2/3, we obtain the results 
given in Table 2. The central value for each p.d.f. is calculated using the median and 
the error corresponds to 34% probability regions on each side of the median. Asymmet- 
ric errors are symmetrized changing the quoted central value. ^ Figures 1 and 2 show, 
respectively, the probability density functions (p.d.f. 's) for some UT parameters and the 
selected region in the p — 77 plane. 



Parameter 


68% 


95% 


99% 


P 


0.196 ± 0.045 


[0.104,0.283] 


[0.073,0.314] 


77 


0.347 ± 0.025 


[0.296,0.396] 


[0.281,0.412] 


«n 


96.1 ± 7.0 


[82.1,110.0] 


[77.7,114.8] 


m 


23.4 ± 1.5 


[20.8,26.1] 


[20.2,27.1] 


in 


60.3 ± 6.8 


[47.0, 74.2] 


[42.5,78.9] 


sin2Q; 


-0.21 ± 0.24 


[-0.65,0.27] 


[-0.77,0.41] 


sin 2/3 


0.726 ± 0.028 


[0.670, 0.780] 


[0.651, 0.797] 


sin(2/3 + 7) 


0.947 ± 0.038 


[0.852, 0.996] 


[0.813,0.998] 


ImAt[10-5] 


13.3 ± 0.9 


[11.5,15.1] 


[10.9,15.6] 



Table 2: Values and probability ranges for the UT parameters obtained from the UT fit 
using the following constraints: \Vuh\ / \Vcb\, Arrid, Arris, £k and sm2p. The value of 
ImAt = ImV^gVtd is also given. 



3.1 Fundamental test of the Standard Model in the quark sector 

The standard fit, illustrated in Figure 2, gives a clear picture of the success of the SM. A 
crucial test consists in establishing CP violation by using the sides of the UT, i.e. CP- 
conserving processes such as the semileptonic B decays and B^^s ~ Bd,s oscillations. The 
comparison of the region selected by these constraints and the one selected by the direct 
measurements of CP violation in the kaon (ek) or in the B (sin 2/3) sectors is shown in 
Figure 3 and gives a picture of the success of the SM in the fiavour sector. This pictorial 

^For the p.d.f. with multiple solutions, e.g. for the UT angles in Section 4, we consider the range 
corresponding to 68% probability and delimited by the intercept of the distribution with an horizontal 
line. 
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Y[°] sin(2p+Y) 

Figure 1: From top to bottom and from left to right, the a-posteriori p.d.f. 's for p, f], a, 
sin 2/3, 7 and sin(2/5 + 7). The red (darker) and the yellow (lighter) zones correspond 
respectively to 68% and 95% of the area. The following constraints have been used in the 
UT fit: \Vub/\Vcb\, ^rrid, ^rris, sk, and sm2j3 from the measurement of the CP asymmetry 
in the J/ipK^ decays. 
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-1 -0.5 0.5 1 _ 

P 

Figure 2: Allowed regions for p and r] using the parameters listed in Table 1. The 
closed contours at 68% and 95% probability are shown. The full lines correspond to 95% 
probability regions for each of the constraints, given respectively by the measurements of 
\Vub\ I \Vch\, Am^, Am<j, Ek, and sin 2/3 from the measurement of the CP asymmetry in 
the J/ipK^ decays. 



agreement is quantified through the comparison between the value of p and rj computed 
from CP-conserving and CP-violating observables: ^ 



p = 0.169 ± 0.057([0.055, 0.310] at 95%) 
1] = 0.364 ± 0.037([0.252, 0.430] at 95%) 



UT sides only 



p = 0.241 ±0.070([0.098, 0.363] at 95%) 1 3.^-^0)+^ U) 
= 0.311±0.030([0.260, 0.371] at 95%) J ^'^ ^ ^ 

Another test can be performed by comparing the value of sin 2/5 from S{J/il)K^) and 
the one determined from "sides" measurements 

sin 2/3 = 0.734 ±0.043 ([0.616,0.811] at 95%) UT sides only 

sin 2(3 = 0.726 ± 0.037 ([0.652, 0.800] at 95%) S{J/tpK^) (5) 

For completeness, we also give the value of sin 2/3 obtained by using all the constraints 
but the direct determination: 

sin 2p = 0.725 ± 0.043 ([0.634, 0.804] at 95%) UT sides + gk- (6) 



^The second allowed region for p and rj from the CP-violating measurements (see Figure 3) has been 
discarded since it is ruled out by the measurement of cos 2/3, see Section 4.3.4. 
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-1 -0.5 0.5 1 _ 

P 

Figure 3: The allowed regions for 'p andrj (close contours at 68%, 95% probability ranges), 
as selected by the measurements of \Vub \ / iKfel? ^"^d? and by the limit on Am^, are com- 
pared with the bands (at 68% and 95% probability ranges) from the measurements of 
CP-violating quantities in the kaon (ek) and in the B (sin2/?j sectors. 



As a matter of fact, the value of sin 2/3 was predicted, before its first direct^ measurement 
was obtained, by using all other available constraints (|V^b| / iKbl, Am^, Am^, and Sk)- 
The "indirect" determination has improved regularly over the years, as shown in Figure 4, 
where the direct measurement is also reported, as a reference. 

The agreement of these determinations confirms the validity of CKM mechanism in the 
SM. This test relies on several non-perturbative techniques, such as the Operator Product 
Expansion for computing B decay rates, the Heavy Quark Effective Theory and LQCD, 
which are used to extract the CKM parameters from the experimental measurements. 
The overall consistency of the UT fit gives confidence on the theoretical tools. Assuming 
the validity of the SM, it is possible to perform a more quantitative test of the non- 
perturbative techniques as discussed in following. 

3.2 Determination of other quantities 

In the previous sections we have obtained the a-posteriori p.d.f.'s for all the UT parame- 
ters. It is also instructive to remove from the fitting procedure the external information 
on the value of one (or more) of the constraints. 

In this section we study the distributions of Am^ and of the hadronic parameters. In 
the first case we do not include in the analysis the experimental information on Bg-Bg 
mixing coming from LEP and SLD. In the case of the hadronic parameters, we remove 

•^In the following, for simplicity, we will denote as "direct" ( "indirect" ) the determination of any given 
quantity from a direct measurement (from the UT fit without using the measurement under considera- 
tion). 
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Figure 4: Evolution of the "indirect" determination o/ sin 2/3 over the years (until 2002). 
From left to right, they correspond to the following papers [4, 8]: dRLMPPPPSOO, 
DDGN90, LMMR92, AL94, CFMRS95, BBL95, AL96, PPRS97, BF97, BPS98, PS98, 
AL99, CFGLM99, CPRS99, M99, CDFLMPRSOO, B.et.al.OO, HLLLOO and CFLPSS [2, 
3]. The dotted lines correspond to the 95% C.L. regions (the only information given in 
those papers). The larger hatched bands (from year '99) correspond to values o/ sin 2/3 
from direct measurements (±lcr). The most recent values contained in this paper are also 
shown. 



from the fit the constraints on their values coming from lattice calculations, and use them 
as free parameters of the fit. In this way we can compare the uncertainty obtained on a 
given quantity through the UT fit to the present theoretical error on the same quantity. 

3.2.1 The expected distribution for Arris 

The p.d.f. for Arris obtained by removing the experimental information coming from Bs- 
Eg mixing is shown in Figure 5. The result of this exercise is given in Table 3 (upper line). 
The present experimental analyses of Eg-Es mixing at LEP and SLD have established a 
sensitivity of 18.3 ps~^ and they show an evidence at about 2(j for a positive signal at 
around 17.5 ps^^, well compatible with the range of the Am^ distribution from the UT fit 
(Figure 5). The inclusion of this information in the UT analysis has a large impact on the 
determination of Am^, as shown by the comparison between the permitted range obtained 
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when this information is either included or not (see Table 3). Accurate measurements of 
Arris, expected from the TeVatron in the near future, will provide an ingredient of the 
utmost importance for testing the SM. 




ArTisEps"'] 

Figure 5: Am^ probability distribution, obtained without using the experimental infor- 
mation from Bg-Bg mixing. 



Parameter 


68% 


95% 


99% 


Am^ (without Am^) [ps^"*^] 


21.2 ± 3.2 


[15.4,27.8] 


[13.8,30.0] 


Am^ (including Am^) [ps^-*^] 


18.5 ± 1.6 


[15.6,23.1] 


[15.1,27.3] 



Table 3: Central values and ranges for Amg corresponding to defined levels of probability, 
obtained by including (or not) the experimental information on Arris . 



3.2.2 Determination of fssy Bbsi Bk and ^ 

To obtain the a-posteriori p.d.f. for a given hadronic quantity, we perform the UT fit 
imposing as input a uniform distribution in a range much larger than the expected interval 
of values assumed for the quantity itself. Table 4 and Figure 6 (left column) show the 
results when one single parameter is taken out of the fit with the above procedure. The 
central value and the error of each of these quantities have to be compared to the current 
evaluations from LQCD, given in Table 1 and also shown in Table 4 for the reader's 
convenience. 

Some conclusions can be drawn. The precision on fs^ \J Bb, obtained from the fit has 
an accuracy which is better than the current evaluation from LQCD. This proves that the 
standard UT fit is, in practice, weakly dependent on the assumed theoretical uncertainty 

on fBsJBB,. 
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Parameter 


68% 


95% 


99% 


C-UTftt 


L15 ± 0.11 


[0.97, 1.44] 


[0.93, 1.59] 


C -LQCD 


1.24 ± 0.06 


[1.13, 1.34] 


[1.10, 1.37] 


/B^y^(MeV) -UTfit 


265 ± 13 


[242, 305] 


[236, 330] 


fn^^jBu.illcV) -LQCD 


276 ± 38 


[201. 350] 


[179. 373] 




0.69 ± 0.10 


[0.53, 0.93] 


[0.49, 1.07] 


Bk -LQCD 


0.86 ± 0.11 


[0.67, 1.05] 


[0.62, 1.10] 



Table 4: Values and probability ranges for the non-perturbative QCD parameters, if the 
external information (input) coming from the theoretical calculation of these parameters 
is not used in the fit. The LQCD values are given with errors obtained by convoluting the 
Gaussian and flat contributions reported separately in Table 1. 



The result on Bk indicates that values of Bk smaller than 0.49 are excluded at 99% 
probability, while large values of Bk are compatible with the prediction of the fit obtained 
using the other constraints. The present estimate of Bk from LQCD, which has a 15% 
relative error (Tabic 1), is as precise as the indirect determination from the UT fit. 

The present best determination of the parameter ^ comes from LQCD. 

In the above exercise we have removed from the fit individual quantities one by one. 
It is also interesting to see what can be obtained taking out two of them simultane- 
ously. Figure 6 shows the regions selected in the planes (/b^ \J Bbs, Bk), {Bk, ^-iid 

(/bs \J Bbs, 0- The corresponding results are summarized in Table 5. For (/s^ \J Bb^, Oj 
the 2-dimensional distribution is not limited by the fit. In such a case, the probability 
attached to a given range depends completely on the ranges of the a-priori distributions 
chosen for the two variables. For this reason, we do not quote any range in Table 5. Still, 

from the plot in Figure 6, a combined lower bound fs^'^ Bb^ ^ 200 MeV and ^ > 0.9 
emerges. 



Parameter 


68% 


95% 


99% 


Bk 


0.68± 0.10 


[0.52, 0.97] 


[0.48, 1.21] 


fB^BBM^y) 


257 ± 15 


[232, 304] 


[221, 348] 


Bk 


0.68 ± 0.18 


[0.45,1.14] 


[0.41,1.33] 




1.31 ± 0.21 


[0.96,1.72] 


[0.89,1.98] 



Table 5: Values and probability ranges for the non-perturbative QCD parameters, if two 
external pieces of information (inputs) coming from the theoretical calculations of these 
parameters are not used in the fit. 
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FBs^[GeV] ^ 



Figure 6: One- and two-dimensional probability distributions for Bk, ^ ond fBs\jBBs 
(see Section 3.2.2). The dark and the light (blue) zones correspond respectively to 68% 
and 95% of the area. The la ranges for these quantities, as calculated from LQCD, are 
indicated with error bars. 
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4 New Constraints from UT angle measurements 



The values for a, 7, sin(2/3+7), and sin 2/3 given in Table 2 have to be taken as predictions 
for future measurements. A strong message is given for instance for the angle 7. Its 
indirect determination is known with an accuracy of about 10%. 

Thanks to the huge statistics collected at the B factories, new CP-violating quantities 
have been recently measured allowing for the direct determination of a, 7, 2/3+7, 
cos 2/3. In the following, we study the constraints induced by these new measurements in 
the p — rj plane and their impact on the global UT fit. 



4.1 Determination of the angle 7 using DK events 



Various methods related to -B ^ DK decays have been proposed to determine the UT 
angle 7 [9]-[ll], using the fact that a charged B can decay into a D^{D )K final state via 
a VcbiYub) mediated process. CP violation occurs if the and the D decay to the same 
final state. These processes are thus sensitive to the phase difference 7 between Vub and 
Vcb- The same argument can be applied to S — > D*K and B — > D^*'>K* decays. 
Three methods have been proposed: 

• Gronau-London-Wyler method (CLW) [9]. It consists in reconstructing the neutral 
D meson in a CP eigenstate: B^ — > D^p^K^, where D^p^ are the CP cigcnstates 
of the D meson. In this case, one can define four quantities, sensitive to the value 
of the angle 7: 

TiB+ D^^p±K+) + T(B- D%p±K-) , 

Rcp± = — ^ = 1 + r| ± 2rBC0S7C0s5B 

V{B+ D^K+) + T{B- D K-) 

^ r{B+ ^ D''cp±K+) - T{B- ^ D''cp±K~) ^ ±2rB sin 7 sin (5^ 

r(B+ ^ D^^p^K+) + T{B- ^ D^^p^K-) Rcp± ^ ^ 

where is the absolute value of the ratio of the Cabibbo-suppressed over the 
Cabibbo-allowed amphtudes: 



rsiDK) 



A{B- D^K- 



The phase 5b in Eq. (7) is the strong phase difference between the h ^ u and the 
b ^ c decays. The main limitation of this method comes from the dilution of the 
interference effect, since tb is expected to be small (~ iKjfc/Kftl)- One can never- 
theless repeat this measurement for several fcp final states of D meson. Moreover, 
this method (and the ones introduced below) can be generalized to D*K, DK* 
and D*K* final states, having different values for and Sb, but the same func- 
tional dependence on 7. In case oi B ^ D*K decays, there is an effective strong- 
interaction phase difference of 180° whenever the D* is reconstructed as D* — > D7 
or D* Dn [12]. 

Atwood-Dunietz-Soni method (ADS) [10]. It consists in forcing the 5° (D^) meson, 
coming from the Cabibbo-suppressed (Cabibbo-allowed) b ^ u {b c) transition to 
decay into the Cabibbo-allowed (Cabibbo-suppressed) Ki: final state. In this way. 
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one can look at the interference between two amplitudes having the same order of 
magnitude. Two quantities are defined: 



R 



ADS — 



^ADS — 



T{B- [K+n-jpK-) + T{B+ ^ [K-tt+]dK+) 
T{B- [K~7r+]DK~) + T{B+ [K+7r-]DK+) 
i"dcs + + ^r-Br-Dcs cos 7 cos{6b + 5d) 

T{B~ [K+tt-]dK-) - T{B+ [K-tt+]dK+) 
T{B- [K+7r-]DK-) + T{B+ [K-7r+]DK+) 
^ ^tbTdcs sin 7 sin(5B + 5d)/Rads (9) 

which are functions (as in the case of the GLW method) of 7, tb and Sb- Since in 
this case and are forced to decay into different final states, the expressions in 
Eq. (9) also depend on the ratio of the Cabibbo-suppressed over the Cabibbo-allowed 
decays of the D, 



and on the additional strong phase shift 5d introduced by these decays. One can 
also use D*K, DK* and D*K*, having different values of and Sb, but the same 
values of 7, vdcs and 6d- 

• Dalitz method [11]. It consists in studying the interference between the h ^ u 
and the h c transitions using the Dalitz plot of D mesons reconstructed into 
three-body final states (such as, for instance, — > Kgii'Ti:^). The advantage 
of this method is that the full sub-resonance structure of the three-body decay is 
considered, including interferences such as those used for GLW and ADS methods 
plus additional interferences due to the overlap between broad resonances in some 
regions of the Dalitz plot. The same analysis is also performed using B D*^K 
decays. The Dalitz analysis has only a twofold discrete ambiguity (7, 7 — tt) and 
not a fourfold ambiguity as in the case of the GLW and ADS methods. One can fit 
the experimental data to extract a 3D likelihood as a function of r^, 7 and 5b- 

Both BaBar and Belle have presented results applying the three methods to several 
final states. It is important to notice that in the case of ADS and Dalitz measurements we 
used directly the experimental fikelihood because of the presence of non-Gaussian effects: 
For Dalitz analyses, there are non-trivial correlations between tb and 7 and the error on 
7 is proportional to l/r^; for ADS analyses, the likelihood of Rads is bounded to be 
positive. 

The p.d.f.'s of 7 and r^, and the selected region in the 7 — plane are shown in 
Figure 7, together with the impact of this measurement on the p — 77 plane. 
The comparison between the direct and the indirect determination is: 

7n = 60.3 ±6.8 ([47.0, 74.2] at 95%) indirect - UTfit (11) 

(59.1 ±16.7 U -120.3 ±17.2 , ^r*w.r*^ 

7 ° = <^ direct -D^ifW 

[([24.7,97.9] U [-155.4, -82.7] at 95%) 

where the two values of the direct determination approximately correspond to the discrete 
ambiguity of the Dalitz method, which dominates the determination of 7. An important 
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-1 -0.5 0.5 



Figure 7: Left column: a-posteriori one-dimensional p. d.f. 's (from top-left to bottom-left) 
for rsiDK), rB{D*K), rB{DK*) and 7, using the results from B — > decays. 
The plots on the right show the measurement of the angle 7 in the 7 — (first three rows) 
and J) — rj planes (bottom). 
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(Jbservable 


DK 


D*K 


DK* 


Acp+ (GLW) 


0.22 ± 0.11 


-0.14 ± 0.18 


-0.07 ± 0.18 


Acp- (GLW) 


0.02 ± 0.12 


0.26 ± 0.26 


-0.16 ± 0.29 


i?cp+(GLW) 


0.91 ± 0.12 


1.25 ± 0.20 


1.77 ± 0.39 


Rcp- (GLW) 


1.02 ± 0.12 


0.94 ± 0.29 


0.76^°J° 


Rads 


0.017 ± 0.009 


< 0.16@90% C.L. 




^ADS 


n 40+0-53 
'-'•^'^-0.46 






rB(Dalitz)-Belle 


0.21 ±0.08 ±0.03 ±0.04 


0.12^°;}? ± 0.02 ± 0.04 




7[°](Dalitz)-Belle 


68 ±15 ±13 ±11 


75±57± 11 ± 11 




r B (Dalitz)-BaBar 


< 0.19 @90% C.L. 


0.155l[J;[J^?± 0.040 ±0.020 




7[°](Dalitz)-BaBar 


70 ± 26 ± 10 ± 10 





Table 6: Summary of the experimental results obtained at the B factories from B 
D^*^K^*^ decays using the GLW, ADS and Dalitz methods [5, 13]. The last two lines give 
the results from the Dalitz analysis, for which the ambiguity 'y tt — 'y is implicit. 




Figure 8: Regions of the plane selected by the constraints imposed by the determination 
of \Vub/Vcb\ and 7 from tree-level processes. 



result of this analysis is also the distributions for the different r^'s: 

rB{DK) = 0.10 ±0.04 ([0.03, 0.17] at 95% C.L.) (12) 
rB{D*K) = 0.09 ± 0.04 ([0.01, 0.18] at 95% C.L.) 
rB{DK*) = 0.32 ± 0.21 ([0.02, 0.74] at 95% C.L.) 
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The small values for rsiDK) and rB{D*K) are limiting the precision on the present 
determination of 7. 

It is interesting to observe that the determination of 7 discussed in this section is 
not affected by new physics (NP) under the mild assumption that NP does not change 
tree-level processes. Therefore, together with the measurement of it already 

provides a constraint on the p-fj plane which must be fulfilled by any NP model. The 
regions selected by these two constraints are shown in Figure 8. We obtain 

p= ±(0.21 ±0.10), ?7 = ±(0.36 ±0.06). (13) 

A more detailed and quantitative analysis of the impact of the UT fit on specific NP 
models will be presented in a forthcoming paper. 

4.2 Determination of 2/^+7 using D^*'>7r{p) events 

The interference effects between the b ^ c and b ^ u decay amplitudes in the time- 
dependent asymmetries of B decaying into 0^*^71 and D^*^ p final states allow for the 
determination of 2/3+7. The time-dependent rates, in the case of Dt: final states, can be 
written as 

(1 + C cos(Amdi) + S sm{Amdt) ) (14) 
(1 — C cos(Amrft) — S sin(Am,;f:) ) 
(1 + C cos{Amdt) - S sm{Amdt) ) 
(1 - C cos(Amdt) +S sm(Amdt) ) 

where S, S and C are defined as 

S = S{D-n+) = 2-^sin(2/? + 7-5) (15) 

1 + 

S = S{D+7r-) = 2^^sin(2/3 + 7 + 5) 

c ^ 







= N 


e 


-rt 


r(b"- 


D-7r+) 




e 


-rt 


R{B^ - 


^ D+n-) 


= N 


e 


-Ft 


R(^- 




= N 


e 


-rt 



1 + 



and r and 5 are the absolute value and the strong phase of the amplitude ratio A{B 
D~7r~^)/ A{B^ — > D~7T^). All these expressions can be generalized to the case of D*7i and 
Dp final states, which have in principle different values of r and S. The ratio r is expected 
to be rather small being of the order of AjT^b/ycbl — 0.02. 

The extraction of the weak phase is made even more difficult by the presence of 
a correlation between the tag side and the reconstruction side in time-dependent CP 
measurements at B factories [19]. This is related to the possibility that the interference 
between b ^ c and b ^ u transitions in i? ^ DX decays occurs also on the tag side. In 
such a case, it is useful to replace S and S by two new parameters a and c which can be 
written as [5] 

S -\-^ 

a = — — ~ 2rsin(2/? + 7)cos((5) (16) 
S — ~S 

c = ~ 2cos(2/5 + 7)(rsin(5) -r'sin(5')) 
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retaining only linear terms in r and r', where r' and 6' are the analogue of r and S 
for the tag side. It is important to stress that the interference in the tag side cannot 
occur when B mesons are tagged using semileptonic decays. In other words, r = when 
only semileptonic decays are used. In the following we will consider the observable a, ci 
(denoting c evaluated for lepton-tagged events), which are functions of r, S and 2/5+7. 
The experimental situation is summarized in Table 7. 



Parameter 


Dn 




Dp 


a 


-0.045 ± 0.027 


-0.030 ± 0.014 


-0.005 ± 0.049 


Q 


-0.035 ± 0.035 


0.010 ± 0.021 


-0.147 ± 0.082 



Table 7: Summary of the experimental results from BaBar and Belle Collaborations. 
World averages are calculated by HFAG [5] from [20]. 




Figure 9: Distribution of r vs. 2/5 + 7 foi" -D*7r decays (top). For illustration, we also 
show the p.d.f. o/2/? + 7 (bottom-left) and the constraint in thep — f] plane (bottom-right) 
obtained assuming SU{3) flavour symmetry and neglecting the annihilation contribution 
m A(B° -> D*-TT+). 

With the present experimental data, a determination of 2/5+7 cannot be obtained from 
D^*^7i(p) modes alone. The number of free parameters exceed the available constraints as 
shown by Eqs.(16). Without further input, one can only find correlations among 2/5+7 
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and the hadronic parameters. For example, the correlation between r{D*7T) and 2/3+7 is 
shown in Figure 9. An independent information on the hadronic parameters would allow 
a determination of 2/3+7. For instance, assuming SU (3) flavour symmetry and neglecting 

annihilation contributions, one can estimate BR{B — > L)*~7r+) from BR{B — > D*'"7r+), 
obtaining r{D*Ti) = 0.015 ± 0.006 ± 0.005 (where the first error is statistical and the 
second is a guessed theoretical error associated to the SU(3) breaking effect and to the 
size of annihilation contributions [20]). Under these two assumptions, we get a constraint 
on 2/3+7 ^ shown in Figure 9. The same can be done for the other two modes. This 
strategy suffers from theoretical uncertainties which cannot be reliably estimated. For 
this reason, we do not use any bound on 2/3+7 the global fit. A more fruitful use of 
these data could be possible when the experimental analyses included the coefficient of the 
cosine term (in such a case the formulae in Eq. (16) should not be expanded in r and r'). 
Additional decay modes, such as S° ^ D^Kg or the Dahtz analysis of B^ — > D^Kgi^^ 
decays [21], would also help. 



4.3 Determination of the angle a 

The angle a can be obtained using the time-dependent analyses of B^ — > tt+tt", B^ — > 
p+p- and 5° {pirf. 

In the absence of contributions from penguin diagrams, these decays give a measure- 
ment of sin 2a. Penguin diagrams with c and t quarks in the loop introduce an additional 
amplitude with a different weak phase. In this case, the experimentally measured quan- 
tity is sin2Q;e//, which is a function of sin2Q; but also of unknown hadronic parameters. 
Several strategies have been proposed to get rid of this so-called "penguin pollution" . 



4.3.1 Isospin analysis oi B ^ tttt and B ^ pp 

Assuming SU(2) flavour symmetry and neglecting electroweak penguins, the decay am- 
plitudes of i3 ^ 7r7r(pp) can be written as 

A^^ = — ^ [e-^"(T + Tc e^-^^O] 
v2 

= — ^ [e-^^Te e^-^^^ + P e'-^^j , (17) 
v2 

where T, and P are real parameters ^ and 5p and 5t^ are strong phases (the strong phase 
of the T term is conventionally set to zero). It should be noted that these parameters 
are different for B ^ Tin and B ^ pp decays. For B pp, following the experimental 
results, only longitudinally polarized final states have been considered. 

In Table 8, we collect the experimental information on these modes. It is important to 
remark the situation of the measurements of S and C in P — > tttt decays. The BaBar and 
Belle Collaborations, even with increased statistics, confirmed the disagreement between 
their results (at about 3.2cr) [16]. For this reason, we exclude these measurements from 

■^These parameters are directly related to the RGI quantities of ref. [23]. In particular, up to a 
trivial rescaling, P is the charming penguin parameter Pi of [23], while T = Ei — pGiM _|_ 
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the global fit waiting for a clarification of the experimental results in the near future and 
use only the experimental informations from B ^ pp decay modes. The result on a from 
pp is shown in Figure 10 and given in Table 9. 





Figure 10: A-posteriori one-dimensional p. d.f.'s of a from pp (top-left), p-K (top-right), 
and their combination (bottom-left) . On bottom right we show the impact of the a deter- 
mination on the p — f] plane. 

The different branching ratios can also be obtained a-posteriori and compared with 
the measured input values. The results of this test are presented in Table 9. In Figure 11 
we show the 2D distribution of the B ^ pp case. 

4.3.2 Dalitz analysis of (p7r)° 

The time-dependent study of B^ (p7r)° decays on the Dalitz plot is a powerful way 
to obtain a [17]. The sensitivity of this analysis depends upon the measurements of the 
branching ratios and CP asymmetries for the various intermediate resonances as well as 



20 







vrvr 






PP 




Observable 


BaBar 


Belle 


Average 


BaBar 


Belle 


Average 


C 


-0.09 ± 0.16 


-0.58 ± 0.17 




-0.23 ± 0.28 




-0.23 ± 0.28 


S 


-0.30 ± 0.17 


-1.00 ± 0.22 




-0.19 ± 0.35 




-0.19 ± 0.35 




4.7 ± 0.6 


4.4 ± 0.7 


4.6 ± 0.4 


30.0 ± 6.0 




30.0 ± 6.0 




5.8 ± 0.7 


5.0 ± 1.3 


5.5 ± 0.6 


22.5 ± 8.1 


31.7 ± 9.8 


26.4 ± 6.4 




1.17 ± 0.34 


2.32 ± 0.53 


1.51 ± 0.28 


0.54 ± 0.41 




0.54 ± 0.41 



Table 8: Experimental inputs from isospin analyses in B ^ tttt and B pp decays [5]. 
The B pp decays are assumed to be fully polarized, in agreement with available mea- 
surements. 



PP 


Output Value 


95% 


99% 


«[1 


96.2 ± 5.5 U 175.3 ±15.4 


[73.7, 122.7] U [153.4, 192.5] 


[68.7,202.7] 


BR{p+p-) 


32.3 ± 5.4 


[21.7,43.0] 


[18.5,46.2] 


5i?(p+pO) 


20.3 ± 4.4 


[11.9,28.8] 


[9.5,31.4] 


BR{p'>p^) 


0.6 ± 0.7 


[0.1,1.4] 


< 1.6 



Table 9: Output values for the main parameters and experimental observable entering the 
isospin analysis. Branching ratios are quoted in units of 10~^. 



DC 



m 



45 ^ 

40 =- 



BR(P°P°) 



brcp^p") 



Figure 11: Plots showing the correlation of BR{B^ —>■ p'^p^) vs BR{B^ p^p^)(left) and 
BR{B^ —>■ p~^p~) vs BR{B^ p^ p^) (right) . The crosses indicate the la experimental 
ranges for the branching ratios, quoted in units of 10~^. 

from the interference between different amplitudes through mixing. In the SM and without 
further theoretical assumptions, one can write the amplitudes of the various intermediate 
decay modes of the B meson as 





= T+-e 


-ia _|_ p-\ — 




= T-+e 


-ia _j_ p — h 


^00 


= TO°e- 


ia _j_ ^00 
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where T and P parameters are complex amplitudes, carrying their own strong phase, 
and the first (second) superscript refers to the p (tt) charge. The amplitudes for the 
CP-conjugate decays can be parameterized in the same way. The two sets of relations 
are characterized by 13 unknowns, one of which is a global phase that can be removed 
and one is fixed by the normalization of the Dalitz plot. Moreover, using SU(2) 
flavour symmetry, can be written as a function of and P^+. This leaves 9 
quantities to be determined including a. Recently, BaBar reported a result of this 
analysis [18] which gives the values of 16 experimental observables that can be written as 
functions of these 9 unknowns. The result for a is shown in the top-right plot of Figure 10. 

4.3.3 Combined results on a 

The results on a from pp, pvr, their combination and the bound on the 'p — rj plane are 
shown in Figure 10. We obtain 

a = (107± 11)° U (176± 14)° ([80, 138] U [150, 197] @95%) (19) 

Notice however that the interpretation of this result can be misleading. Indeed, if an 
independent information would allow discarding the second solution at 176°, the error 
associated to the first determination of a would be larger than 11°. 

4.3.4 Determination of cos 2/3 from JtjjK*^ decays 

From the time-dependent analysis of the decay J/ipK*'^, it is possible to extract 

both sin 2/3 and cos 2/3 [27]. The results obtained by Belle and BaBar are barely compatible 
(see Table 10) and we combined them using the skeptical approach of ref. [28], assuming 
the values S = 1.3 and A = 0.6 for the parameters of the model. We also varied these two 
parameters without obtaining sizable deviations for the observed combined p.d.f. 

cos 2/3 

BaBar 3.32l°:^^±0.27 
Belle 0.31 ± 0.91 ± 0.11 
Skeptical 1.9 ± 1.3 



Table 10: Experimental values for cos 2/3 and our skeptical average. 

The skeptical likelihood of cos 2/3 and the impact of this determination on the p — fj 
plane, once the a-priori bound |cos2/3| < 1 is imposed, are shown in Figure 12. The 
probability of cos 2/3 being less than zero is about 13%. This implies that this measurement 
removes the ambiguity associated to sin 2/3, suppressing one of the two allowed bands for 
sin 2/3 in Figure 2. 

4.4 Determination of the Unitarity Triangle parameters using 
also the new UT angle measurements 

It is interesting to sec the selected region in p — r/ plane from the direct measurements of 
the UT angles: sin 2/3, 7, a, and cos 2/3. The plot is shown in Figure 13. In Table 11 we 
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Figure 12: The skeptical likelihood o/ cos 2/3 (left) and the cos 2/5 constraint on the p — rj 
plane (right). 




P 



Figure 13: Allowed regions for J) andrj obtained using the measurements of the UT angles 
only: sin 2/5, a, 7, and cos 2/5. The closed contours at 68% and 95% probability are shown. 
The full zones correspond to 95% probability regions from individual constraints. 



report the results obtained using these constraints. 

The results given in Table 12 are obtained using all the available constraints: 
\Vub\ I iKfel, ^f^di Am^, sin 2/5, 7, a, and cos 2/5. Figure 14 shows the corresponding 
selected region in the J) — f] plane. 

Given the present experimental uncertainties, the new measurements of UT angles, 
taken individually, would loosely constrain the values of p and f]. In addition, some depen- 
dence of the results on the choice of the a-priori distributions is present, particularly for 
those quantities which are poorly constrained by the experiments (as it should be in the 
Bayesian approach). On the other hand, when combined with the sin 2/5 measurement. 
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Figure 14: Allowed regions for p andfj using the parameters listed in Table 1 together with 
the new UT angle measurements. The closed contours at 68% and 95% probability are 
shown. The full lines correspond to 95% probability regions for each of the constraints, 
given by the measurements of \Vub\ / \Vcb\, Arrid, Am^, Ek, sin2/5, 7, a, and cos2/?, re- 
spectively. 



they select an area in the 'p-fj plane comparable to the one available in the pre-S-factory 
era and largely independent of the chosen a-priori distributions. The agreement between 
this area and the output of the standard UT fit is an important demonstration of the 
consistency of the CKM mechanism in describing non-leptonic B decays and CP asym- 
metries. 

5 Compatibility plots 

In this section we discuss the interest of measuring with a better precision the various 
physical quantities entering the UT analysis. We investigate, in particular, to which extent 
future and improved determinations of the experimental constraints, such as sin 2/?, Arris 
and 7, could allow us to possibly invalidate the SM, thus signalling the presence of NP 
effects. 

5.1 Compatibility between individual constraints. The pull dis- 
tributions. 

In CKM fits based on a minimization, a conventional evaluation of compatibility 
stems automatically from the value of the at its minimum. The compatibility between 
constraints in the Bayesian approach is simply done by comparing two different p.d.f.'s. 

Let us consider, for instance, two p.d.f.'s for a given quantity obtained from the UT 
fit, /(xi), and from a direct measurement, /(X2): their compatibility is evaluated by 
constructing the p.d.f. of the difference variable, X2 — xi, and by estimating the distance 
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Parameter 


68% 


95% 


99% 


P 


0.241 ± 0.081 


[0.090,0.441] 


[0.036,0.753] 


V 


0.328 ± 0.038 


[0.249,0.411] 


[0.182,0.468] 


«n 


102 ± 11 


[78, 124] 


[66, 134] 


pn 


23.6 ± 1.8 


[19.9,26.9] 


[18.2,28.7] 


in 


54.1 ± 11.7 


[30.6,76.6] 


[21.9,85.6] 


sin2Q; 


-0.38 ± 0.36 


[-0.94,0.38] 


[-0.99,0.66] 


sin 2/3 


0.727 ± 0.037 


[0.654,0.800] 


[0.632,0.823] 


sm{2p + 7) 


0.955 ± 0.038 


[0.805,0.998] 


[0.159,1.0] 


ImAt[10-5 


12.8 ± 1.6 


[9.6, 16.0] 


[8.2, 17.3] 



Table 11: Values and probability ranges for the UT parameters obtained by using the 
constraints: sin 2/3, 7, a, and cos 2/5. 



Parameter 


68% 


95% 


99% 


P 


0.207 ± 0.038 


[0.129,0.282] 


[0.106,0.308] 


V 


0.341 ± 0.023 


[0.296, 0.386] 


[0.282, 0.400] 


an 


97.9 ± 6.0 


[86.0, 109.7] 


[82.5,113.7] 


m 


23.4 ± 1.5 


[20.7,26.1] 


[20.2,27.1] 


7[°] 


58.5 ± 5.8 


[47.3, 70.2] 


[43.5,73.8] 


sin2Q; 


-0.27 ± 0.20 


[-0.64,0.13] 


[-0.75, 0.25] 


sin 2/3 


0.725 ± 0.028 


[0.669,0.779] 


[0.651,0.795] 


sm{2p + 7) 


0.958 ± 0.030 


[0.884,0.997] 


[0.853,0.998] 


ImAtflO-s; 


13.2 ± 0.9 


[11.5,14.8] 


[11.0,15.3] 



Tabic 12: Values and probability ranges for the Unitarity Triangle parameters obtained by 
using all the available constraints: \Vub\ / \Vcb\, Am^, ArUg, Sr, sin 2/3, 7, a, and cos 2/3. 

of the most probable value from zero in units of standard deviations. The latter is done by 
integrating this p.d.f. between zero and the most probable value and converting it into the 
equivalent number of standard deviations for a gaussian distribution ^ . The advantage of 
this approach is that no approximation is made on the shape of p.d.f. 's. In the following 
analysis, f{xi) is the p.d.f. predicted by the UT fit while the p.d.f of the measured 
quantity, f{x2), is taken Gaussian for simplicity. The number of standard deviations 
between the measured value, X2 ± a{x2), and the predicted value (distributed according 
to f{xij) is plotted as a function of X2 (x-axis) and a{x2) (y-axis). The compatibility 
between xi and X2 can be then directly estimated on the plot, for any central value and 
error of the measurement of X2- 

^In the case of Gaussian distributions for both xi and X2; this quantity coincides with the pull, which 
is defined as the difference between the central values of the two distributions divided by the sum in 
quadrature of the r.m.s of the distributions themselves. 
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If two constraints turn out to be incompatible, further investigation is necessary to 
tell if this originates from a "wrong" evaluation of the input parameters or from a NP 
contribution. 

5.2 Pull distribution for sin 2/3 

We start this analysis by considering the measurement of sin2/3. The left plot in Figure 15 
shows the compatibility ("pull") between the measurement of sin 2/5 and its indirect de- 
termination, obtained in the SM using the constraints |Kib/|K;fe|, Am^, Am^, and Ek (but 
excluding S{J/ijjK^)), as function of the measured value (x-axis) and error (y-axis) of 
sin 2/?. The cross indicates the present experimental average of S{J /ipK*^). 




sin2|3 sin2|3 



Figure 15: The compatibility ("pull") between the direct and indirect determination of 
sin 2/? as a function of the measured values and errors. On the left, the indirect dis- 
tribution of sin 2 P is computed from the UT fit including \Vub/Vcb\, Am^, Am^, and Sk, 
without using the direct measurement from — >■ J/ipK^ . On the right, the indirect distri- 
bution includes also the constraint from the measurements o/ sin 2/5 from J/ipK^. 
The compatibility regions from la to 6a are displayed. The crosses display the posi- 
tion (value/error) of the measurements o/ sin 2/5 from B^ —>■ J /ipK^ (left plot) and from 
B — *■ (f)K^ (right plot), extracted from the HFAG average for S{(f)K^). We also mark 
with an asterisk the value of sin 2(3 obtained using instead our skeptical combination for 

The plot shows that, considering the present precision of 0.037 on the measured value 
of sin 2/?, the 3a compatibility region is in the range [0.51, 0.88]. Values outside this range 
would be, therefore, not compatible with the SM prediction at more than 3a level. To get 
these values, however, the presently measured central value should shift by at least 6a. 

The conclusion that can be derived from Figure 15 is the following: although the 
improvement on the error on sm2(3 has an important impact on the accuracy of the UT 
parameter determination, it is very unlikely that in the near future we will be able to use 
this measurement to detect any failure of the SM, unless the other constraints entering the 
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fit improve substantially or, of course, in case the central value of the direct measurement 
moves away from the present one by several standard deviations. 



5.2.1 Role of sin 2/3 from Penguin processes 

It was pointed out some time ago that the comparison of the time-dependent CP asymme- 
tries in various B decay modes could provide evidence of NP in B decay amplitudes [22] . 
Since sin 2f3 is known from S{J /ipK^), a significant deviation of the time-dependent asym- 
metry parameters of penguin dominated channels from their expected values would indi- 
cate the presence of NP. 

Prom this point of view, the pure penguin h — > sss processes as S — > are the 
cleanest probes of NP. The only SM uncertainty in the extraction of sin 2/3 from S{(j)K'^) 
comes from the penguin matrix elements of current-current operators containing up-type 
quarks. This is expected to be small and can be estimated in a given model of hadron 
dynamics and constrained using the experimental data. For example, using the model of 
ref. [25], we obtain 



In the previous equation we have used both the HFAG average for S{(f)K^) and our es- 
timate obtained using the skeptical approach of ref. [28] (with S = 0.6 and A = 1.3). 
The average value of S{(f)K^) is obtained using four measurements having /NDF=2.Q. 
As this result is used to find evidence for NP it has been considered that systematic 
uncertainties may have been underestimated, so errors have been inflated using the ap- 
proach [28]. If the present discrepancy were instead due to a statistical fluctuation, it 
should disappear in the future and the two approaches will converge to the same result. 

The compatibility of sin 2/9 shown in the right plot of Figure 15 has been obtained 
using all the constraints of the standard analysis, namely iKtft/Kbl, ^^d, Am^, Sk, in- 
cluding sin 2/3. The cross and star on this plot correspond to the values of sin 2/3 in 
Eq. (20), extracted from the HFAG and the skeptical average respectively. Depending 
on the chosen average procedure, the compatibility of sin 2/3 from B — > (f)K^ with the 
indirect determination ranges from about la to 2(t. 

The extraction of sin 2/3 could be extended to channels receiving additional contribu- 
tions from b — > uus transitions, such as S — > rj'K^, B^ — > foK^ or S — > n'^K'^. However, 
in these cases, hadronic uncertainties are more difficult to estimate and expected to be 
channel-dependent. For this reason, CP asymmetries in these channels can be different 
and should not be naively averaged. A detailed analysis of these modes goes beyond the 
goal of this work and will be presented elsewhere. 

5.3 Pull distribution for Anis 

The plots in Figure 16 show the compatibility of the indirect determination of Am^ with 
a future determination of the same quantity, obtained using or ignoring the experimental 
information coming from the present bound. 

From the plot in Figure 16 we conclude that, once a measurement of Arris with an 
expected accuracy of ~ 1 ps~^ is available, a value of Am^ greater than 32 ps~^ would 
imply NP at 3cr level or more. 




0.34 ± 0.20 ± 0.08 from S^ko = 0.34 ± 0.20 (HFAG) 
0.46 ± 0.28 ± 0.08 from S^k^ = 0.46 ± 0.28 (skeptical) 



(20) 
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Figure 16: The compatibility between the direct and indirect determination of Anis, as 
a function of the value of Ams(ps^^), using (left) or ignoring (right) the present experi- 
mental bound. 

5.4 Pull distribution for the angles 7 and a 




y[°] an 

Figure 17: The compatibility between the direct and indirect determination of'-/ (left) and 
a (right), as a function of their measured value and corresponding uncertainty using the 
UTfit results. The cross displays the position (value/error) of the present measurements 
of 1 from charged B mesons decaying into D*K* final states. 

The left plot in Figure 17 shows the compatibility of the direct and indirect deter- 
mination of 7. It can be noted that, even in case the angle 7 can be measured with a 
precision of 10° from B decays, the predicted So" region is still rather large, corresponding 



28 



to the interval [25,/, 95]°. Values larger than 100° would clearly indicate physics beyond 
the Standard Model. 

The present direct determination of the angle 7, being in perfect agreement with UT 
fit, cannot provide any evidence of NP independently of the precision of the measurement. 
To a lesser extent, the same conclusion can be drawn for a, see right plot in Figure 17. 
Assuming again a precision of 10°, NP at 3a shows up outside the range [58,/, 132]°. 

6 Conclusions 

Flavour physics in the quark sector has entered its mature age. Today the Unitarity 
Triangle parameters are known with good precision. A crucial test has been already 
done: the comparison between the Unitarity Triangle parameters, as determined with 
quantities sensitive to the sides of the triangle (semileptonic B decays and oscillations), 
and the measurements of CP violation in the kaon (e^-) and in the B (sin2/9) sectors. The 
agreement is "unfortunately" excellent. The Standard Model is "Standardissimo" : it is 
also working in the flavour sector. This test of the SM has been allowed by the impressive 
improvements achieved on non-perturbative methods which have been used to extract the 
CKM parameters. 

Many B decay branching ratios and CP asymmetries have been measured at B facto- 
ries. The outstanding result is the determination of sin 2/3 from B hadronic decays into 
charmonium-X^ final states. On the other hand many other exclusive hadronic rare B 
decays have been measured and constitute a gold mine for weak and hadronic physics, 
allowing in principle to extract different combinations of the Unitarity Triangle angles. 

Besides presenting an update of the standard Ulfit analysis, we have shown in this 
paper that new measurements at B factories begin to have an impact on the overall picture 
of the Unitarity Triangle determination. In particular the angle 7 is today measured 
through charged B decays into DK final states within 20% accuracy and only a twofold 
ambiguity. In the following years the precise measurements of the UT angles will provide 
further tests of the Standard Model in the flavour sector to an accuracy up to the percent 
level. 

Finally, by introducing the compatibility plots, we have studied the impact of future 
measurements for testing the SM and looking for new physics. In the near future the 
measurement of Am^ and of the UT angle 7 will play the leading role. 
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